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A single Yb atom is loaded into a high-finesse optical cavity with a moving lattice, and its nu- 
clear spin state is manipulated using a nuclear magnetic resonance technique. A highly reliable 
quantum state control with fidelity and purity greater than 0.98 and 0.96, respectively, is confirmed 
by the full quantum state tomography; a projective measurement with high speed (500 /is) and 
high efficiency (0.98) is accomplished using the cavity QED technique. Because a hyperfine cou- 
pling is induced only when the projective measurement is operational, the long coherence times 
(Ti = 0.49 s and T2 = 0.10 s) are maintained. Our technique can be applied for implementing a 
scalable one-way quantum computation with a cluster state in an optical lattice. 

PACS numbers: 03.65.Wj, 03.67.Mn, 32.80.-t, 42.50.Dv 



Oh: 



in 
> 

00 

rn 

o 
o 



X 



Quantum computation enables us to solve NP prob- 
lems such as factorization into prime factors and to per- 
form quantum simulations of solid-state models'll]. Al- 
though various research activities have been dedicated 
to the implementation of quantum computing, actual ex- 
periments have so far been limited to only few steps of 
quantum operations using a small number of qubits. Re- 
cently, an array of quantum gates was implemented for 
neutral atoms with massively parallel operation using 
optical lattice systems [3-i4j. However, in these experi- 
ments, the coherence time was limited to the order of 
milliseconds, which is too short to achieve scalable quan- 
tum computation[5|, |6[. The other challenge is to address 
single sites of the lattice and control single qubits in a 
reliable manner. 

It is advantageous to implement a qubit with a sin- 
gle nuclear spin|7Hll| because the nuclear spin is robust 
against a stray magnetic field; this robustness can be 
attributed to the weak strength of the magnetic mo- 
ment generated by the nuclear spin as compared with 
that generated by the electronic spin 12]. Here, we em- 
ploy a single nuclear spin of the ^''^Yb atom in the form 
of a qubit in an optical lattice potential and implement 
both manipulation and state verification of the qubit with 
high fidelity and purity using nuclear magnetic resonance 
(NMR) and cavity QED techniques. It should be noted 
that a long coherence time of the qubit originates from 
an extremely weak coupling between the nuclear spin and 
the electromagnetic field and is normally traded off for 
enabling fast and highly efficient state measurement. We 
couple the electronic spin with the nuclear spin using 
hyperfine interaction in an excited state, and such a cou- 
pling is induced only when we irradiate the atoms with 
the probe light used for readout. We can thus maintain 
a long coherence time of the qubit. In the cluster com- 



puting, both high readout efficiency and high-speed de- 
tection are mandatory; these two conditions are satisfied 
in our experiment by using a cavity QED setupflSi, ll^J. 
We will discuss the cluster quantum computation using 
the cavity QED system in the last section. 

Our experimental setup, relevant energy levels, and 
the time chart of the experiment are shown in Fig. 1. 
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FIG. 1: Experimental setup and time chart, (a) Energy-level 
diagram of ^^^Yb. (b) Experiment setup, (c) Time chart for 
the experiment. |t) and |4,) denote the magnetic sublevels 
ruF = +1/2 and — 1/2 in the ground state ^So(F = 1/2), re- 
spectively. The substates in the excited state "^Pi are labeled 
as Ittt-f). 



We construct a qubit with a 1/2 nuclear spin of the 
ground state "^So(F = 1/2) of the ^'^^Yb atom and de- 
fine the quantized axis along the cavity axis by using the 
bias magnetic field induced by "Coill." In the following 
discussion, |t) and \\.) represent the magnetic sublevels of 
^So(F = 1/2). Details of Zeeman slowing and the double 
MOT system have been obtained from Refs. 17 and 18. 
To load the atoms into the cavity mode, we utilize the 
concept of a "moving lattice." We first activate an opti- 
cal standing wave potential, which is created by counter- 
propagating frequency-doubled YVO4 laser beams (532 
nm), with a peak potential depth of 110 fiK . Because the 
optical standing wave overlaps the MOT beams, atoms 
are loaded to the lattice potential. After the MOT beams 
and the relevant magnetic field are turned off, the lattice 
potential starts moving with a speed that is proportional 
to the value of S. By changing the value of S accord- 
ing to the relation 6(t) = Sq sin(7ri/T) , atoms are trans- 
ported to the cavity mode, where Sq = 2tt x 700 kHz and 
r = 100 ms . To initialize the nuclear spin state to |^) , 
we irradiate the atoms in the moving lattice with a cir- 
cularly polarized spin polarization beam that is resonant 
with the iSo(F = 1/2) ^ iPi(F' = 1/2) transition (399 
nm). When we manipulate the nuclear spin state, the 
radio frequency (RF) of magnetic field is applied to the 
atoms by using "Coil2" after the atoms arrive in the cav- 
ity mode, where the radio frequency is resonant with the 
Zeeman splitting between two substates |t) and \l) (Fig. 
la, 6g ~ 2tt X 2.5 kHz ). Because of the RF irradiation, 
the nuclear spins in the atoms rotate around the axis or- 
thogonal to the cavity axis. Our cavity consists of two 
concave mirrors having ultra-high reflectivity at 556 nm 
and is characterized by the following three parameters: 
the maximum interaction rate between atoms and pho- 
tons go = 2tt X 2.8 MHz , cavity decay rate (HWHM of 
the cavity resonance line) k — 2tt x 4.8 MHz , and atom 
decay rate (half natural linewidth) 7 = 27r x 91 kHz . 
The length of our cavity is stabilized to 150 /im by in- 
jecting a 560-nm laser beam and utilizing the FM side- 
band method such that the cavity is resonant with the 
iSo(F = 1/2) -^ 3p^(-p/ ^ 3/2) transition (556 nm). It 
should be noted that the far-off resonant locking beam 
(560 nm) is responsible for trapping the atoms in the 
cavity mode, where the trap depth is 30 /iK (beam waist 
is 19 /im ). To perform the projective measurement, we 
increase the bias magnetic field and the Zeeman split- 
ting in the excited state '^Pi(F' = 3/2), where the split- 
ting 6e ~ 2tt X 60 MHz is much larger than the cavity- 
enhanced linewidth F = 7 [l + 2.9^/(^7)] = 27rx3.4 MHz 
for the iSo(F = 1/2) -^ ^p^^p/ ^ 3/2) transition. When 
the atoms are irradiated with a linearly polarized pump 
beam resonant with the iSo(F = 1/2) -)■ 3Pi(F' = 3/2) 
transition, only the \D atoms are excited in a cyclic man- 
ner and repeatedly emit fiuorescence photons into the 
cavity mode. Therefore, the detection of more than one 
photon implies that the nuclear spin is projected to H) 
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FIG. 2: Rabi oscillation and Ramsey interference. Time de- 
velopments of population of \\) state as measured by Rabi 
oscillation (a) and Ramsey interference (b). The error bars 
are estimated on the basis of the statistical distribution of 
signal counts. 



. To perforin experiments with single atoms, we prepare 
dilute atoms such that the expected mean atom num- 
ber in the cavity mode becomes less than unity; utilize 
the post-selection method to perform the experiments. 
The total detection efficiency of a photon emitted from 
an atom is 0.1; this value is enhanced by the Purcell ef- 
fect. We obtain four photons from each single atom on 
average. On the basis of the assumption of the binomial 
distribution for photodetection, the efficiency of the pro- 
jective measurement is estimated to be 0.98. The result 
obtained from the projective measurement contain errors 
that can be attributed to the dark counts of the detector 
and also to unwanted spin fiips caused by the excitation 
of the It) state. These lead to about 2% error for the 
diagonal elements of the density matrix. 

The time chart (Fig. Ic) of the experiment consists 
of two parts, i.e., the preparation and detection of the 
nuclear spin states. Arbitrary single nuclear spin states 
are prepared by rotating nuclear spins around two or- 
thogonal axes. The rotation around the cavity axis is 
accomplished using the Larmor precession generated by 
the bias magnetic field, and the other rotation is imple- 
mented using Rabi oscillation through irradiation of the 
resonant RF field. Figure 2a shows the Rabi oscillations 
of the single nuclear spin generated by the RF field. We 
induce Zeeman splitting between the two ground states 
It) and It) using the bias magnetic field and apply a 
square-shaped of the RF pulse whose frequency is reso- 
nant with the splitting S — 2Trx2.5 kHz. The population 
of the down spin state as a function of the RF irradia- 
tion time is shown in Fig. 2a. Experimental results show 
that the visibility of 0.96, which is estimated by sinu- 
soidal fitting. This high visibility enables us to precisely 
control the nuclear spin rotation around one axis. The 
fact that the operation is carried out using single nuclear 
spins can be confirmed as follows. In our projective mea- 
surement, we consider that the spin is projected to |t) 
when the photodetector counts more than one photon. 
Suppose that there are two atoms in the cavity mode. In 
this case, the photon counts are obtained even if the spin 
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FIG. 3: Full quantum state tomography. Real and imagi- 
nary parts of density matrices for the prepared nuclear spin 
states reconstructed using the maximum likelihood estima- 
tion. The corresponding spin states are (a) (|t) + |i))/-\/2 , 
(b) (lt)+ili))/V2,and(c) !;). 



state is |||) , ||t) ; or |t|) , and thus, the variation in 
the probability of detecting more than one photon is no 
longer sinusoidal (see a dotted curve in Fig. 2a). The 
curve obtained from the experimental results represents 
a simple sinusoidal curve, which indicates the successful 
post-selection of a single nuclear spin. 

To check the controllability of two-axis rotations, we 
measure the Ramsey interference for a single nuclear spin. 
On the basis of the measurements of the Rabi oscillations 
(Fig. 2a), we use an RF field with a duration of 3.2 ms as 
a 7r/2 pulse. We irradiate single nuclear spins with two 
7r/2 RF pulses, and the nuclear spins rotate around the 
cavity axis during the interval between the two pulses. 
The population of the down spin state as a function of 
the time interval between the two RF pulses is shown 
in Fig. 2b. The curve obtained from the experimental 
results is in good agreement with the theoretical curve, 
and the visibility is calculated from the sinusoidal fitting 
to be 0.99, which is consistent with the value expected 
from the error generated in the projective measurement. 
It should be noted that the obtained visibility of Ramsey 
interference (0.99) is slightly higher than that obtained 
from the Rabi oscillations. In our experiment, the sta- 
bility of the bias magnetic field is better than that of the 
RF field. We believe that the above difference arises be- 
cause the RF irradiation time in the case of the Ramsey 
interference experiment is shorter than that in the case of 
the Rabi oscillation experiment. Hereafter, we use only 
7r/2 pulses of the RF field. By appropriately varying the 
timing of the 7r/2 pulses, superposition states and also 
eigenstates are prepared. 

By using spin polarization and rotation, we prepare 
arbitrary states of single nuclear spins. Here three typi- 
cal spin states are prepared, i.e., (a) (|t) + |i)) /v^ , (b) 
(It) + * ID) /\/2 , and (c) ||) ; these sates are defined on a 
rotating frame with a frequency of Sg determined by the 
bias magnetic field. At t = ms, we begin transporting 



single nuclear spins from the MOT to the cavity with the 
moving lattice. At t = 50 ms, the nuclear spin is polar- 
ized to the down state using the spin polarization beam 
(399 nm). Superposition states are created by applying 
the 7r/2 pulse at (a) t = 100.0 ms and (b) t = 100.1 ms. 
To reconstruct the density matrix (p) for each state, we 
perform projective measurements for the single nuclear 
spins along three directions. To perform the projective 
measurements along three orthogonal axes, we rotate the 
spins using Rabi oscillation or Lamor precession tech- 
niques before irradiating the atoms with the pump beam. 
Figure 3 shows the density matrix (p) for each state re- 
constructed by using the maximum likelihood estimation 
method. The purity {p = Tr [p^] ) and the fidelity to the 
ideal state (f — {ipi\p |V'i); IV'i) is the ideal state) are es- 
timated to be (a) (p,f) = (0.98 ± 0.01, 0.99 ± 0.005) , 
(b) (p,f) = (0.96 ±0.01, 0.98 ±0.005) and (c) (p,f) = 
(0.97 ± 0.03, 0.98 ± 0.02) ; the Monte Carlo method is 
employed for this estimation assuming a binomial distri- 
bution for the signal counts. The purities and fidelities 
for the Fig. 3 are slightly less than 1, which are limited by 
the distortion of the RF field and photon scattering due 
to the moving lattice. One-way quantum computing re- 
quires a highly efficient projective measurement because 
the computation is based on the result of individual pro- 
jection and the scalability is therefore determined by the 
detection efficiency. Efficient and fast detection are real- 
ized in our experiment by using hyperfine interaction and 
also the enhanced mode-selective spontaneous emission 
due to the Purcell effect derived from the cavity QED 
system [ig|. Although the detection efficiency is deter- 
mined by the rate of photon count (4 counts/500 /is) in 
our experiment, a dramatic improvement (0.9998) can be 
expected by simply collecting atomic fluorescence from 
both sides of the cavity. Such an extremely high value of 
efficiency enables us to perform thousands of operations 
in the one-way quantum computing. 

To evaluate the longitudinal relaxation time Ti, we 
first measure the lifetime r of single Yb atoms in the 
micro-cavity mode. We prepare the spin state of \l.) and 
measure the population of |J,) as a function of the trap- 
ping time. Next, the spin state of |t) is prepared, and 
the dependence of the population of H) on the trapping 
time is measured. The sum of them gives the survival 
probability of single Yb atoms in the micro-cavity mode 
as a function of the trapping time (see the inset of Fig. 
4a); the hfetime r is estimated to be 0.44 ± 0.03 s by 
fitting the data to an exponentially decaying function. 
Figure 4a shows the population of \\.) as a function of the 
trapping time; the fitting is obtained by normalizing the 
population of \D with the survival probability. From the 
fitting, the longitudinal relaxation time Ti is estimated 
to be 0.49 ±0.15 s. To estimate the transverse relaxation 
time T2, we measure the visibility of the Ramsey inter- 
ference as a function of the trapping time (Fig. 4b). The 
insets of Fig. 4b show the typical Ramsey interferences, 
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FIG. 4: Measurement of coherence time of single nuclear 
qubit. (a) Decay of population of |^) state. The inset shows 
the survival probability of the atom in the cavity as a function 
of time, (b) Time development of visibility of Ramsey fringe. 
The inset shows typical examples of Ramsey interference pat- 
terns. For both cases (a) and (b), the error bars are evaluated 
on the basis of the binomial distribution of signal counts. 



where the horizontal axis represents the time separation 
between two RF pulses and the vertical axis represents 
the population of the down spin state. The transverse 
relaxation time T2 is estimated to be 0.10 ± 0.01 s from 
the fitting. It should be noted that in the case of neutral 
alkali atoms, the coherence time of Zeeman sublevels is 
only of the order of 10 /xs[7|. Even in the case in which 
both the clock states of different hyperfine substates and 
the spin echo technique are used, the coherence time is 
limited to 30 ms[8|. The coherence time of 0.1 s obtained 
in our experiment is thus much longer than that of alkali 
atoms. Here, we discuss the possible sources responsible 
for the decoherence. The experimentally obtained values 
of Ti and T2 are limited by photon scattering due to the 
moving lattice and by the fluctuation or inhomogeneity 
of the bias magnetic field. Because the fluctuation of 
the bias magnetic field does not affect the population 
of the spin state, Ti is considered to be determined by 
the photon scattering rate Tp according to the relation 
Ti = l/Fp . The photon scattering rate Tp is estimated 
to be 2.0 s~^ from the experimental parameters including 
the imperfection of mode match of counter propagating 
beams. The inverse of Tp is 0.5 s which is in agreement 
with the experimentally obtained value of Ti = 0.49 s. 



Here, we define the decoherence rate Tm , which origi- 
nates from the energy-conserving dephasing effect due to 
the fluctuation and inhomogeneity of the bias magnetic 
field. Because transverse relaxation occurs due to longi- 
tudinal relaxation and the dephasing effect, the relation 
I/T2 = 1/Ti + r„j should be satisfied. From the experi- 
mentally obtained values of Ti and T2, F^ is estimated 
to be 8 s~^, which corresponds to the magnetic field fluc- 
tuation of 9 mG. The ratio of the coherence time and the 
time required for the projective measurement is 200; this 
ratio approximately gives the possible operation number. 
The suppression of the magnetic field fluctuation below 
1 mG is technically possible; suppression to such values 
will increase the coherence time to values up to 0.5 s. 

Our cavity QED system can be applied to implement 
a large scale of the cluster quantum computation. We 
prepare the nuclear spin qubit per site of an optical lat- 
tice using the fermionic band insulator state. To create 
the macroscopic cluster state for the nuclear spins, the 
nuclear spin dependent potential is required. While such 
a potential can be created by a laser beam whose de- 
tuning is less than hyperfine splitting, the ratio of the 
natural linewidth to the hyperfine splitting is not enough 
to suppress the spontaneous emission in the presence of 
ordinary dipolar allowed transitions of alkaline atoms. 
However, the problem can be solved by making an opti- 
cal potential with the ^Sq — > '^P2 ultra-narrow transition 
of ytterbium. To create the cluster state, we use the 
Ising-type interaction based on the magnetic dipole mo- 
ment of '^P2 state because the scattering length of the 
^^^Yb atom is too small. Since the ■^P2 state is auto- 
matically mixed when we induce the optical potential, 
the required interaction can be realized by simply mak- 
ing close the distance between neighboring atoms by the 
nuclear spin dependent potential. The required phase 
fiip is estimated to take 10 ms. Individual addressing 
is implemented with a magnetic gradient and a specific 
atom is excited to the ^P2 state using the Rabi oscilla- 
tion. The nuclear spin state is transferred to the spin 
state of the Zeeman sub-levels of '^P2 state. The atom in 
the '^P2 state is transported to the micro cavity with our 
belt-conveyer technique and it is returned to the nuclear 
spin state in the micro cavity. We can perform the selec- 
tion of basis and projective measurement on this qubit 
by using the method presented in this article. The clus- 
ter computation is thus accomplished by repeating these 
processes. 
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